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Band gap engineering in finite elongated graphene nanoribbon heterojunctions: 

Tight-binding model 

Benjamin O. Tayo 

Physics Department, Pittsburg State University, Pittsburg, KS 66762, USA 

A simple model based on the divide and conquer rule and tight-binding (TB) approximation is 
employed for studying the role of finite size effect on the electronic properties of elongated graphene 
nanoribbon (GNR) heterojunctions. In our model, the GNR heterojunction is divided into three 
parts: a left (L) part, middle (M) part, and right (R) part. The left part is a GNR of width Wl, 
the middle part is a GNR of width Wm, and the right part is a GNR of width Wr. We assume that 
the left and right parts of the GNR heterojunction interact with the middle part only. Under this 
approximation, the Hamiltonian of the system can be expressed as a block tridiagonal matrix. The 
matrix elements of the tridiagonal matrix are computed using real space nearest neighbor orthogonal 
TB approximation. The electronic structure of the GNR heterojunction is analyzed by computing 
the density of states. We demonstrate that for heterojunctions for which Wl = Wr, the band gap 
of the system can be tuned continuously by varying the length of the middle part, thus providing a 
new approach to band gap engineering in GNRs. Our TB results were compared with calculations 
employing divide and conquer rule in combination with density functional theory (DFT) and were 
found to agree nicely. 


I. INTRODUCTION 

Graphene is a two-dimensional (2D) allotrope of car¬ 
bon with excellent electronic and mechanical properties, 
making it suitable for multiple applications in nanoscale 
electronics and nanophotonicsA major deficiency in 
graphene’s properties is the absence of a band gap ren¬ 
dering it impossible for use in switching circuitsi^ Sev¬ 
eral approaches have been used to induce a band gap in 
graphene such as electrically gated bilayer graphene 
substrate induced band gapjii^ or isoelectronic codoping 
with boron and nitrogen^ Recently, it has become pos¬ 
sible to engineer the band gap of graphene by etching 
or patterning along a given direction to produce ultra 
narrow quasi one-dimensional (ID) nano sheets referred 
to as graphene nanoribbons (GNRs)i2r— with remark¬ 
able properties such as width-dependent tunable band 
gaps, exciton-dominated optical spectra, and room tem¬ 
perature ballistic transport As quasi ID materials, 
GNRs are extremely sensitive to their surrounding con¬ 
ditions, which provides a route for manipulating their 
electronic properties. Additionally, other factors such as 
finite size effectedge effectfi^^— and the presence of 
strain^i^— could be used to effectively tune the electronic 
properties GNRs. 

As exciting as the properties of GNRs are, a ma¬ 
jor bottleneck in exploiting their properties is the dif¬ 
ficulty of synthesizing high quality GNRs having spe¬ 
cific widths and edges. The most common method 
for synthesizing GNRs is lithographic patterning out of 
larger 2D graphene sheets. Such a top-down approach 
lacks precision and often results in GNRs with random 
shapes, rough edges and poor electronic quality. Re¬ 
cently, a group in Berkeley Lab and UG Berkeley have 
successfully demonstrated a precise method for synthe¬ 
sizing high quality GNRs with controlled structure via 
bottom-up molecular engineering^^. They successfully 


synthesized ultra-thin graphene nanoribbon heterojunc¬ 
tions with non-uniform width, providing a new approach 
to band gap engineering. Motivated by such studies and 
the potential of molecular bottom-up engineering in pro¬ 
ducing ultrathin GNR heterojunctions with controlled 
widths and finite lengths, we decided to carry out de¬ 
tailed studies on the electronic properties of finite length 
ultrathin GNR heterojunctions using TB approximation. 
TB approximation is a very simple, reliable, and compu¬ 
tationally efficient method for electronic structure calcu¬ 
lations and has been used successfully for modelling the 
electronic, optical and transport properties of carbon- 
based systems 

In this paper, we present a model based on the divide 
and conquer rule^^ and TB approximation for studying 
the role of finite size effect on the electronic properties 
of elongated GNR heterojunctions. In our model, the 
GNR heterojunction is divided into a left (L) part, mid¬ 
dle (M) part, and right (R) part. The left part is a GNR 
of width Wl, the middle part is a GNR of width Wm, 
and the right part is a GNR of width Wr. We assume 
that the left and right parts of the GNR heterojunction 
interacts with the middle part only so that there is no in¬ 
teraction between the left part and the right part. Under 
this approximation, the total Hamiltonian of the system 
can be expressed as a block tridiagonal matrix. The indi¬ 
vidual matrix elements of the total Hamiltonian operator 
describing the left, middle and right parts, as well as the 
coupling matrix between right and middle and left and 
middle can readily be computed using the real space near¬ 
est neighbor orthogonal TB approximation. The elec¬ 
tronic structure of the GNR heterojunction is analyzed 
by computing the density of states. We demonstrate that 
for heterojunctions for which Wl = Wr, the band gap 
of the system can be tuned continuously by changing the 
length of the middle part, thus providing useful informa¬ 
tion for experimental developments. Because our model 
is based on the TB approximation, significant compu- 
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type atoms and N B-type atoms, as shown in Fig. [T] 
The dangling cr-bonds at the edges are passivated by H 
atoms (or other atoms/groups like O and OH). Conse¬ 
quently, the C-C bond length at the edges is shortened 
to (1 — Sc)ac, where Sc (3 to 5%) is the compressive strain 
on the edge C-C bond due to H passivation^i^i^. Ad¬ 
ditionally, N-AGNRs can be classified into three distinct 
families JV = 3p, 3p-|-l, 3p-|-2, where p is a positive integer 
and their electronic properties are known to exhibit dis¬ 
tinct family splittingi^^— For the finite extended AGNR 
(Fig. nib)), we introduce an integer length parameter 
such that the length L of the system is L = N^T. We will 
define an N-M GNR heterojunction as a heterojunction 
formed between an N-AGNR and an M-AGNR. 


A. Infinite Periodic GNRs 


FIG. 1. (a) Infinite periodic H-passivated AGNR showing the 
number of dimer lines along the width of the ribbon, (b) 
Schematic representation of the finite extended GNR, con¬ 
structed by a replication of the repeating unit (m). m rep¬ 
resents the unit cell of an N-AGNR. (c) Two-leg ladder with 
N rings representing the equivalent TB Hamiltonian of the 
periodic system at the F point. t| and are model pa¬ 

rameters (hopping integrals) that are used in constructing the 
TB Hamiltonian matrix. 


tational efficiency is gained without compromising accu¬ 
racy. 

This paper is organized as follows: In Sec. im we de¬ 
scribe the general formalism. First, we discuss the TB 
approximation for the infinite periodic system (Sec. HIl 
A), then we describe how the divide and conquer rule 
can be applied to describe the properties of finite elon¬ 
gated systems (Sec. HlB), and in Sec. m C we describe 
the GNR heterojunction using a three-part model and 
the block tridiagonal matrix. In Sec. uni we discuss the 
model parameters used and describe how the TB real 
space matrices can be constructed within the divide and 
conquer rule. In Sec. El we discuss the results. A short 
summary concludes the paper. 


II. GENERAL FORMALISM 

We consider an infinite periodic armchair GNR 
(AGNR) of width W = — l)oc and translation 

period T = 3ac, where N is the number of dimer lines 
and Qc 1.43 A the unstrained carbon to carbon (G-G) 
bond length at the center of the GNR (see Fig. UKa)). 
Since the width of an AGNR is specified by the number 
of dimer lines along the ribbon, we will use the notation 
N-AGNR to refer to an AGNR with N dimer lines along 
the ribbon. The unit cell of an N-AGNR contains N A- 


The electronic states of infinite periodic (p) GNRs are 
expressed in terms of the axial momentum (k) and the 
lateral momentum (fc„), where n in an integer describing 
the quantization of the component of electron’s momen¬ 
tum along the width of the ribbon. AGNRs are semicon¬ 
ductors with a direct band gap at the F point. At fc = 0, 
the TB Hamiltonian for an AGNR reduces to a two-leg 
ladder lattice systen:>i^, as shown in Fig. [T] (c). The 
Hamiltonian of this simpler model Tip (where p stands 
for periodic) reduces to 
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where (q/i) denote a site, site energies, 

t\ the nearest neighbor hopping integrals within each leg 
and between the legs respectively, and the annihila¬ 
tion operator of 7r-electrons on the i-th site of the /r-th 
leg. We remark here that in this model, the electronic 
properties of GNRs are sensitive only to the three pa¬ 
rameters: the site energies and the nearest neigh¬ 
bor hopping integrals 4- These parameters 

will differ for perfectly terminating, edge passivated and 
strained GNRs. In general for k ^Q^TLp can be expressed 
in matrix form for the translational invariant system. If 
we order the basis as Ai, B 2 , A 3 , ..., A^r-i, Bj^, and 
Bi, A 2 , B 3 , ..., Bn-i, Aat, then the nearest neighbor 
Hamiltonian can be split into four N x N blocks 
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Here, dk = with T being the lattice constant. 

The electronic band structure of the AGNR can then be 
obtained by solving the eigenvalue equation 


7ip{k)Cxn{k) = Exn{k)Cxn{k), (4) 

where A = c (v) corresponds to the conduction (valence) 
band, and n the band index. The coefficients Cxn{k) are 
TB wave function amplitudes. 


B. Finite Extended GNRs 

The eigenstates of a finite extended GNR are discrete 
and can be labelled by the quantum number i. Finite 
size effects in GNRs are studied using the divide and 
conquer rule. In this approach, the finite-length GNR is 
constructed by a replication of a repeating unit (to), as 
shown in Fig. [T] (b). The TB Hamiltonian T-Lj (where 
f stands for finite size) of the system is then given by a 
block-tridiagonal matrix: 



( Hm 

Hm,m 

0 




H* 

m,m 

Hm 

Hm,m 

0 
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where Hm is the TB Hamiltonian of the subunit to, and 
Hm,m the coupling between nearest neighbor subunits. 
The simplicity of our approach lies in the fact that both 
Hm and Hm,m are approximated using the nearest neigh¬ 
bor TB approximation. Our TB description makes cal¬ 
culations easy to compute compared to other methods 
that use density functional theory (DFT) calculations to 
compute the matrices Hm and Hm,rrr^^- 




7-13 GNR heterojunction 



FIG. 2. (a) Bottom-up synthesis of graphene nanoribbons 

from molecular building blocks 1 and 2. The resulting rib¬ 
bon, or heterojunction, has varied widths as a result of dif¬ 
ferent width molecules 1 and 2. Adapted by permission from 
Macmillan Publishers Ltd: Nature Nanotechnology (Ref. 28), 
copyright 2015. (b) Schematic representation of the equiva¬ 
lent system consisting of three parts: left (L), middle (M), 
and right (R). In this case, the left and right parts are in¬ 
finitely long 7-AGNR, while the middle part is a finite length 
13-AGNR. 


assume that the GNR heterojunction can be divided into 
three parts: the left (L) part which is an N-AGNR, the 
middle (M) part which is an M-AGNR and the right (R) 
part which is an N-AGNR. Again using the divide and 
conquer rule, the TB Hamiltonian Hh (where h stands 
for heterojunction) of the GNR heterojunction can be 
described by the tridiagonal matrix: 

(HfiL) Ulm 0 \ 

dih — ( Hm,l HfiM) 'Hm,r I ■ ( 6 ) 

\ 0 nR,M HfiR) ) 

where 'Hf{L) is the finite length TB Hamiltonian matrix 
of an N-AGNR, 'Hf{M) is the finite length TB Hamil¬ 
tonian matrix of an M-AGNR, and T-Ll,m the coupling 
matrix between the N-AGNR and M-AGNR forming the 
heterojunction. 


C. GNR Heterojunctions 


III. MODEL PARAMETERS AND 
COMPUTATIONAL DETAILS 


We consider a finite elongated GNR heterojunction 
formed between GNRs of different widths (Fig. [2|). We 


In this section, we provide some details about the dif¬ 
ferent parameters used in our model and describe how to 
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calculate the block matrices using the TB approximation. 
For perfectly terminating AGNR, we will set the nearest 
neighbor C-C TB hopping integral to t = 2.7 eV, a value 
that has been used to successfully describe the electronic 
properties of graphene^. For H-passivated AGNRs, the 
bond lengths parallel to dimer lines at edges are com¬ 
pressed by about 3.5% relative to those in the middle of 
the ribbon. A decrease in G-G bond length will increase 
overlap of tt orbitals which leads to an increase in the hop¬ 
ping integral. Likewise, an increase in bond length will 
result to a decrease in tt orbital overlap, which accord¬ 
ingly decreases the hopping integral. The analytic ex¬ 
pressions for TB matrix elements between carbon atoms 
as a function of the C-C bond length can be expressed in 
terms of the Chebyshev polynomials Tm{x) yielding^ 

„ r — 

= ^ CmTm-iiy) -y = (7) 

m=l 2 


where r G (a, b) is the interatomic distance for C-C in¬ 
teractions, and (a, b) the range of values over which the 
expansion is valid. The coefficients and boundaries 
a and b are tabulated in Ref. 41. Using Eq. © , we 
can show that a 3.5% compressive strain on the bond 
length at the edges induces a 12 % increase (5*) in the 
hopping integral . The effect of H-passivation can then 
be accounted for by setting t■ = (1 + St)t, for i = 1 and 
i = N, t for i = 2,..., TV — 1, and = t for 

fi = 1,2, i = l,...iV — 1. Previous studies carried out 
for edge passivated GNRs have shown that to first-order, 
the change in onsite energy due to edge passivation does 
not alter the band gapj ^^’^^ We will therefore assume 
that changes in onsite energies due to H-passivation are 
negligible. Hence we shall set all the onsite energies at 
£fi,i = 0 for /r = 1, 2 and n = 1, 2,..., iV. 

Using the approximations given above, we can show 
that for the finite system, the orthogonal TB nearest 
neighbor (NN) real space Hamiltonian Hm for each sub¬ 
unit is given by 

{ t, if j are NN interior bonds 

(1 -|- 6 t)t, if i,j are NN edge bonds ( 8 ) 
0 , otherwise 

where St — 12 % is the percent increase of the hoping 
integral due to 3.5% decrease in C-C bond length at the 
edges. Using the same arguments and assuming that the 
coupling between replicating units only involve interior 
bonds, the real space TB matrix Hm,m describing the 
coupling between adjacent m units is given by 


m.m 


t, a i,j are NN interior bonds 
0 , otherwise 


(9) 


For all the GNRs considered, we will assume the C-C 
bond length at the center is Oc = 1.43 A and the C-C 
bond length at the edges is (1 — Sc)ac, with Sc = 3.5%. 
For electronic band structure calculations, we will ne¬ 
glect the change in band structure of the AGNR due 
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FIG. 3. DOS of a H-passivated 7-AGNR plotted for different 
values of the ribbon length L = NlT, where T = 0.43 nm 
is the length of the repeating unit used in constructing the 
finite length system. The DOS is computed with an artificial 
broadening of 10 meV and compared to the DOS of the infinite 
periodic system (red dashed curve superimposed on the black 
curve in each panel). The Fermi energy is set to Ef — 0. 


to quasiparticle effects. Quasiparticle effects are known 
to provide significant improvements in the band gap '^^d^ . 
Our TB treatment however provides excellent qualitative 
results that capture the properties of GNRs. Quantita¬ 
tively speaking, the TB approximation yields results that 
are comparable to results obtained using DFT in the local 
density approximation (LDA)^. The advantage of using 
TB model lies in the computational simplicity, compared 
to DFT methods which are very expensive for describ¬ 
ing finite elongated extended systems. All computations 
were carried out using Mathematica software^. 


IV. RESULTS AND DISCUSSION 

We consider a 7-AGNR and an 11-AGNR as choice of 
materials to be studied. The 7-AGNR belongs to the 
3p -h 1 family and has a large band gap, while the 11- 
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FIG. 4. DOS of a H-passivated 11-AGNR plotted for different 
values of the ribbon length L — NlT, where T = 0.43 nm 
is the length of the repeating unit used in constructing the 
finite length system. The DOS is computed with an artificial 
broadening of 10 meV and compared to the DOS of the infinite 
periodic system (red dashed curve superimposed on the black 
curve in each panel). The Fermi energy is set to FIf = 0. 


FIG. 5. DOS of a H-passivated 7-11 GNR heterojunction 
plotted for different values of the length L = NlT of the 11- 
AGNR connecting the left and right segments. The DOS is 
computed with an artificial broadening of 10 meV. The Fermi 
energy is set to Ep = 0. As the length of the 11-AGNR is 
varied, the band gap can be tuned continuously between the 
band gaps of the 11-AGNR and the 7-AGNR, i.e., between 
0.16 eV and 1.54 eV. 


AGNR belongs to the 3p-|-2 family and has a smaller band 
gap. Our choice of material selection is thus reasonable 
since our goal is to study band gap engineering in GNR 
heterojunctions. The 7-11 GNR heterojunction will thus 
provide a significant range of band gaps lying between 
the smaller band gap of the 11-AGNR and the larger 
band gap of the 7-AGNR. Since our goal is to study the 
properties of finite elongated GNR heterojunctions, we 
will start by using our TB formalism to study the role of 
finite size on the electronic properties of the 7-AGNR and 
the 11-AGNR. Then we will study the effect of finite size 
on the 7-11 GNR hetero junction and relate its properties 
to those of the 7-AGNR and an 11-AGNR from which it 
is made. 

We will study the influence of finite size effect on elec¬ 
tronic band structure by plotting the density of states 
(DOS). The finite temperature DOS per electron of the 


infinite periodic (p) system Pp{E) is given by^i^ : 

1 ^ r^/T 

PF(i^) = ^EE/ dk5[E-EUk)im 

where E\n{k) are eigenenergies corresponding to Bloch 
functions, Ne is the total number of tt electrons in the 
GNR, n = 2tt/NlT is the length of the ID reciprocal 
space for each allowed state, Np is the length parameter 
describing the number of repeating units in the N-AGNR 
of finite length, L = NpT, T being the unit cell length. 
The density of states of the finite (f) extended system 
Pf{E) is given by 

1 

PfiE) = J^^Y.^^E-Si), ( 11 ) 

i—1 

where Si are eigenenergies of the discrete system. For 
comparison purposes, both Pp{E) and Pf{E) are normal- 
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FIG. 6. Band gap of a H-passivated 7-11 GNR heterojunction 
plotted for different values of the length L — NlT of the 11- 
AGNR connecting the left and right segments. Dashed black 
lines represent the band gaps of the 7-AGNR (1.54 eV) and 
that of the 11-AGNR (0.16 eV). As the length of the 11- 
AGNR is increased, the band gap decreases monotonically 
from 1.54 eV to 0.16 eV. 

ized such that 

/ OO nOO 

Pp{E)dE= / pf{E)dE = l. 

-OO j —OO 

For computational purposes, we replace the Dirac delta 
function with a Lorentzian with line width F = 0.01 eV. 
We present the DOS for an energy range of ±2 eV around 
the Fermi energy Ep = 0 for several ribbon lengths L — 
N^T, where T = 0.43 nm is the unit cell length. 

Figure [3] and [4] show the DOS as a function of ribbon 
length for a hydrogen passivated 7-AGNR and 11-AGNR, 
respectively. It can be seen that as the length of the 
system increases, the DOS evolves from a set of irregu¬ 
larly spaced energy levels (typical of finite length molec¬ 
ular systems) to van Hove singularities (VHSs) which is 
an important characteristic of ID infinitely periodic sys¬ 
tems. When the length of the ribbon exceeds 101 T, 
finite size effects become negligible and the DOS of the 
finite system converges to the DOS of the infinite peri¬ 
odic system. Similar results have been obtained using 
divide and conquer method in combination with DFT— . 
The advantage of our method lies in the simplicity of the 
model and computational efficiency. The band gap of 
the infinite system can be obtained from the DOS. For 
the 7-AGNR, the band gap is 1.54 eV while that for the 
11-GNR is 0.16 eV. These values agree very well with 
DFT-LDA band gaps^. 

Next we consider the 7-11 GNR heterojunction. The 
left part of the heterojunction is a 7-AGNR of length L = 
101 r, the middle part is an 11-AGNR of variable length 
L = NpT, and the right part is a 7-AGNR of length L — 
101 T. This model allows us to study the properties of the 
7-11 GNR heterojunction as a function of the width of the 
11-AGNR segment, connecting the left and right parts of 
the GNR heterojunction. Figure [S] shows the DOS of the 


hydrogen passivated 7-11 GNR heterojunction plotted as 
a function of the length of the 11-AGNR which forms the 
middle of the heterojunction. Three important regimes 
can clearly be seen from the DOS plot. When Np < 7, 
the DOS is dominated by the DOS of the 7-AGNR and 
the band gap of the system is approximately equal to 
1.54 eV, which is the band gap of the 7-AGNR. When 
7 < Nl < 101, the DOS contains features from both the 
7-AGNR and the 11-AGNR with the band gap varying 
between the band gap of the 11-AGNR and the 7-AGNR, 
i.e., between 0.16 and 1.54 eV. When Np > 101, the DOS 
of the 7-11 GNR heterojunction is dominated by the DOS 
of the 11-AGNR and the band gap converges to 0.16 eV. 
This is a remarkable result and it clearly shows that band 
gap engineering in GNR can be achieved by synthesizing 
GNR heterojunctions with varied widths. 

In Fig. [BJ we plot the band gap of the 7-11 GNR het¬ 
erojunction for different values of the length L = NpT 
of the 11-AGNR connecting the left and right segments. 
The figure shows that by increasing the length of the 11- 
AGNR, we can tune the bang gap of the material con¬ 
tinuously between the band gap of the 7-AGNR (1.54 
eV) and the band gap of the 11-AGNR (0.16 eV). When 
Nl < 7, the band gap of the 7-11 GNR heterojunction 
is approximately equal to that of the 7-AGNR (1.54 eV). 
When 7 < Np <101, the band gap decreases monotoni¬ 
cally from 1.54 eV to 0.16 eV. In this region, the energy 
gap {Eg) varies with Nl as Eg = A + BN£^, where 
A = 0.06 eV and B — 10.6 eV. The fit is good to within 
a mean absolute error of 0.008 eV. When Np > 101, the 
band gap of the system converges to the band gap of an 
11-AGNR (0.16 eV). 


V. CONCLUSION 

In summary, we have shown that the relevance of fi¬ 
nite size effects on the electronic properties of GNRs can 
well be accounted for by using the divide and conquer 
rule combined with the TB approximation. Because our 
model is based on the TB approximation, significant com¬ 
putational efficiency is gained without compromising ac¬ 
curacy. By computing and analyzing the DOS, we found 
that for H-passivated AGNRs, the transition from a dis¬ 
crete molecular system to a continuum periodic system 
occurs at a length of about 43 nm, which is quite con¬ 
sistent with previous results obtained using divide and 
conquer method in combination with DFT—. For the 
7-11 GNR heterojunction, the DOS and energy band 
gap plots as a function of the length parameter Np of 
the middle part reveal three important regimes: when 
Nl < 7, the band gap of the 7-11 GNR heterojunction 
is approximately equal to that of the 7-AGNR (1.54 eV), 
when 7 < Nl < 101, the band gap decreases monotoni¬ 
cally from 1.54 eV to 0.16 eV, and when > 101, the 
band gap of the system converges to the band gap of an 
11-AGNR (0.16 eV). Our studies clearly show that band 
gap engineering in GNRs can be achieved by synthesizing 
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GNR heterojunctions with varied widths, thus providing 
useful information for experimental developments. 
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